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Archean cratonic mantle recycled at a mid-ocean ridge 
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Basalts and mantle peridotites of mid-ocean ridges are thought to sample Earth’s upper mantle. Osmium isotopes 
of abyssal peridotites uniquely preserve melt extraction events throughout Earth history, but existing records 
only indicate ages up to ~2 billion years (Ga) ago. Thus, the memory of the suspected large volumes of mantle 
lithosphere that existed in Archean time (>2.5 Ga) has apparently been lost somehow. We report abyssal peridotites 
with melt-depletion ages up to 2.8 Ga, documented by extremely unradiogenic '8’0s/1*°0Os ratios (to as low as 
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0.1095) and refractory major elements that compositionally resemble the deep keels of Archean cratons. These 
oceanic rocks were thus derived from the once-extensive Archean continental keels that have been dislodged and 
recycled back into the mantle, the feasibility of which we confirm with numerical modeling. This unexpected 
connection between young oceanic and ancient continental lithosphere indicates an underappreciated degree of 


compositional recycling over time. 


INTRODUCTION 
The asthenosphere, a mechanically weak layer below the drifting 
rigid lithospheric plates, is important for mantle convection and 
crust formation on Earth. Although the asthenosphere has been 
widely thought to be well homogenized by mantle convection (1), 
there is now abundant evidence that it is highly heterogeneous in 
composition, containing widely distributed, severely melt-depleted 
mantle domains (2-5). Nevertheless, such ultradepleted geochemical 
signatures are only rarely detected in mid-ocean ridge basalts 
(MORBs) (6) and ocean island basalts (7). This is likely because an 
ultradepleted mantle is highly refractory due to ancient melt 
extraction (5), and it can produce only very small volumes of magma 
that are masked by the large volume of normal MORB generated in 
the melting column. Ancient ultradepleted signals are therefore best 
preserved in residual abyssal peridotites, where their depletion ages 
can be recognized by Os isotopes (8) and Hf isotopes (9). Abyssal 
peridotites with Paleoproterozoic Os model ages as old as 2.2 billion 
years (Ga) have been found (3, 5), but the origin of such ancient 
mantle domains remains controversial. They could be an ancient 
oceanic lithospheric mantle recycled through subduction (5), or 
they might be fragments of a refractory subcontinental lithospheric 
mantle (SCLM) (10). Disruption of the SCLM beneath Archean 
cratons is not uncommon (11, 12), and ancient lithospheric material 
is therefore expected to exist in the asthenosphere. However, abyssal 
peridotites with Archean depletion ages (>2.5 Ga) have been absent 
from the published record. Moreover, if disrupted SCLM can be 
transported by the asthenosphere to mid-ocean ridges forming new 
oceanic lithosphere that is ultimately subducted, the existence of 
ancient SCLM in the convecting mantle at large may be greatly 
underestimated at present. 

The ultraslow-spreading Southwest Indian Ridge (SWIR), with a 
near-constant full spreading rate of ~14 mm year, is an ideal place 
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to test for the potential existence of Archean mantle at a mid-ocean 
ridge. It is located at a distance of ~2000 km from the nearest 
Archean craton, Kaapvaal, in South Africa (Fig. 1A). The SWIR 
spans 7700 km between the Bouvet and Rodriguez triple junctions 
with an overall unusually steep trend (~60°) to the spreading direc- 
tion relative to the ridge axis. It generally has a crustal thickness 
of ~2 to 6 km, with large, widely spaced volcanic edifices and inter- 
vening regions where the mantle is exposed directly on the seafloor 
(13). Named for the cruises on which the samples were dredged, the 
Knorr and Protea peridotites were sampled west and east of the 
Andrew Bain fracture zone (FZ), respectively (Fig. 1, B and C). 
Details of the locations, petrography, and major element geochemistry 
of the studied samples are provided in Materials and Methods. 


RESULTS 

Major and trace elements 

The studied SWIR abyssal peridotites with protogranular to 
porphyroclastic textures have variable compositions (table S1). 
They are all partially serpentinized and weathered, with variable 
values of loss on ignition (LOI) up to 13 weight % (wt %). A few 
Protea samples are quite fresh, with LOI values of <2 wt %. The 
Knorr samples overall have higher SiO) contents than the Protea 
samples (fig. $2). The Protea samples have relatively fertile compo- 
sitions with 1.04 to 3.33 wt % Al.O3 (Fig. 2A), similar to abyssal 
peridotites globally (14). The Knorr samples have more variable 
Al,O3 contents of 0.34 to 2.59 wt % and can be subdivided into 
high- and low-Fo groups (data $2) according to olivine forsterite 
[Fo; 100 x Mg/(Mg + Fe)] contents. The low-Fo Knorr group has 
relatively fertile compositions (ie., with Al,O3 contents of >1 wt %) 
similar to the Protea samples, whereas the high-Fo Knorr (HFK) 
samples are more refractory, with mostly Al,O3 contents of <1 wt %. 
While both the Protea and low-Fo Knorr peridotites have olivine 
Fo of <91.5, the HFK peridotites have Fo of >91.5 (Fig. 2B); the 
latter also have higher spinel Cr# [Cr/(Cr + Al)] values than the 
former (fig. $2). Compared to orthopyroxene (Opx) in the Protea 
and low-Fo Knorr peridotites, Opx in the HFK samples have higher 
Mg# [Mg/(Mg + Fe)] but lower Al,O3 contents (fig. $3, A and B). 
Overall, the HFK samples are the most refractory peridotites ever 
recorded from the seafloor. 
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Fig. 1. Map of the SWIR and sample locations. (A) Spreading directions and half rates are given in millimeters on the arrows. The red dashed arrows are the mantle flow 
field constrained by seismic anisotropy (36). (B) The super oblique segment of the SWIR (9° E to 16°E). The dredge haul KN162-9-37, from which both high-Fo and low-Fo 
peridotites were sampled, is located to the right of the largest non-hot spot volcano, Joseph Mayes Seamount. (C) The central segment of the SWIR (28°E to 42°E), where 
dredge hauls PROT5S-19 and PROT5-24 are located. Other dredged samples shown are from the PetDB database (www.earthchem.org/petdb). TJ, triple junction. 
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Fig. 2. Geochemistry of SWIR abyssal peridotites. (A) Whole-rock MgO versus Al203. (B) Olivine Fo contents versus spinel Cr# values. Fields of cratonic and 
off-craton peridotites are from (56). Data sources of abyssal peridotites, forearc peridotites, and South African (SA) cratonic peridotites are from the PetDB database 


(www.earthchem.org/petdb). 


Both the Protea and low-Fo Knorr samples are relatively fertile; 
their low olivine Fo contents and spinel Cr# indicate a relatively low 
degree of partial melting. Using the function [F = 10 x In(Cr#) + 24] 
between spinel Cr# values and melting starting from a least-depleted 
abyssal peridotite (15), the Protea samples have experienced a 
minimum of 4 to 14% degrees of fractional melting, and the low-Fo 
Knorr samples have experienced 8 to 14% (table S2). Similar extents 
are also inferred from the heavy rare earth element (HREE) con- 
tents of their clinopyroxene (Cpx) (fig. $4 and table S3). The 
relatively fertile compositions of these samples, based on spinel Cr#, 
are not distinguishable from those of global abyssal peridotites (fig. $2). 
By contrast, the HFK samples are more refractory with much higher 
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spinel Cr# values (0.42 to 0.63) and olivine Fo contents (91.5 to 
92.6); their Opx and Cpx both have higher Mg# values (table $2), 
and their estimated degree of fractional melting is a minimum of 
15 to 20%. Such high extents of melting are also indicated by 
remarkably lower Ti and HREE contents in their Opx than those of 
both the Protea and low-Fo Knorr peridotites (fig. S5 and table S3). 
Cpx large enough for laser ablation analysis was absent in almost all 
HFK peridotites except KN162-9-37-66, which is more enriched in 
light rare earth elements (LREEs) relative to HREE (fig. $4), with a 
(Ce/Yb)y [chondrite-normalized according to (16)] ratio of ~8. 
Such a LREE-enriched pattern is atypical of abyssal peridotites (14) 
but common in metasomatized cratonic mantle xenoliths (17). 


2o0f9 


ZZOT ‘OT Jaquiaydag uo Z10°a0uNTOS' MMM //:sdyYy Wo papeojyumog 


SCIENCE ADVANCES | RESEARCH ARTICLE 


Highly siderophile elements and Re-Os isotopes 

Both highly siderophile element (HSE; including Os, Ir, Ru, Pt, Pd, 
and Re) contents and Re-Os isotopes are listed in table $4. The Protea 
and low-Fo Knorr samples have relatively flat HSE patterns (fig. S6A) 
similar to global abyssal peridotites (8, 18), whereas the HFK 
samples are strongly depleted in Pd (fig. S6B). Most HFK samples 
have Pd/Ir ratios remarkably lower than abyssal peridotites globally 
(fig. S7) but similar to cratonic mantle xenoliths that underwent 
high extents of melting (19). The Protea samples have highly variable 
Re contents [i.e., 0.07 to 1.42 parts per billion (ppb)], with some 
higher than the inferred value of the primitive upper mantle [PUM; 
(20); 0.35 + 0.06 ppb]. Their '*’Os/'*"Os ratios range from 0.1163 to 
0.13054, yielding Re depletion ages (Tgp) relative to PUM younger 
than 1.89 Ga. The low-Fo Knorr samples have relatively radiogenic 
'87Qs/'88Os ratios of 0.12288 to 0.13849; most have '*”Os/'**Os 
ratios higher than the PUM, possibly due to seafloor weathering 
(21) or melt refertilization (22), and cannot yield meaningful Os 
model ages. Two of them, with '*’Os/'**Os ratios lower than PUM, 
yield Trp ages younger than 1 Ga. Similarly, young Trp ages (0.23 to 
0.69 Ga) have been reported for five peridotites from the same 
oblique ridge segment (8). In comparison, all HFK peridotites 
except three (KN162-9-37-46, -63, and -84) have very unradiogenic 
'87Qs/'88Qs ratios of 0.1095 to 0.11323, which yield much older Trp 
ages, varying from Paleoproterozoic to Neoarchean (2.28 to 2.78 Ga). 
The HSE patterns of the Knorr samples also indicate the late addi- 
tion of Re (fig. S6B), resulting in their variable Re contents. Never- 
theless, the Knorr peridotites, particularly the HFK peridotites, 
have low '®’Re/'®8Os ratios (commonly <0.1), which have negligible 
effects on their '*’Os/'**Os ratios. Neither the Knorr nor the Protea 
samples exhibit a correlation between whole-rock Al,O3 and 
'87Q5/'88Qs ratios (fig. S8). 


DISCUSSION 

Overall, the HFK peridotites have highly refractory compositions 
that are different from most abyssal peridotites. They have high 
olivine Fo and spinel Cr#; in particular, two samples (KN162-9-37-15 
and -28) have spinel Cr# values higher than 0.6, which are exceedingly 
rare in residual abyssal peridotites (14). High-Cr# spinel occurs 
both in forearc peridotites that experienced hydrous melting (23) 
and in Archean cratonic mantle xenoliths that have experienced 
very high degrees of melting (17). The HFK samples have olivine 
grains (Fo of >91.5) that are more magnesian than those in forearc 
peridotites (23, 24), plotting within the cratonic mantle xenolith 
field (Fig. 2B). It should be noted that the HFK peridotites have 
relatively higher whole-rock SiO, contents than the low-Fo Knorr and 
Protea peridotites (fig. $1), which might reflect the Si-metasomatism, 
a phenomenon that has been previously described for the Archean 
cratonic mantle xenoliths (25). 

The Re-Os isotope system is the best tool to constrain the timing 
of melt extraction experienced by mantle peridotites, and the rhenium 
depletion age (Tgp) provides the minimum estimate of their ages 
(see Materials and Methods). The Protea peridotites show a histo- 
gram of '*’Os/'**Os ratios very similar to abyssal peridotites globally 
(Fig. 3), with Trp ages younger than 1.9 Ga; similarly, old Os model 
ages have been reported for abyssal peridotites from the Mid-Atlantic 
Ridge (2, 3) and the Arctic Gakkel Ridge (5). In contrast, the Knorr 
peridotites show a histogram of '*’Os/!*8Os ratios similar to 
mantle xenoliths from South African cratons. While the low-Fo 
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Fig. 3. Whole-rock '®”0s/"®8Os ratios of SWIR abyssal peridotites. Two dashed 
lines refer to Trp ages of 2.5 and 1.8 Ga, which were calculated relative to the PUM 
(20). Data of abyssal peridotites are from (2-5, 8,27), and data of mantle xenoliths 
from Kaapvaal and Zimbabwe cratons are from (19, 57-64). The higher 18795/188Qs5 
ratios of SWIR peridotites than those of the PUM might be due to interaction with 
seawater (27). 


Knorr peridotites also have relatively young Trp ages not older than 
1 Ga, the majority (17 of 20) of HFK peridotites have notably 
lower '*’Os/'**Os ratios (Fig. 3), corresponding to much older, 
Paleoproterozoic to Archean Trp ages (2.28 to 2.78 Ga). Such Archean 
Trp ages have not been previously reported from abyssal peridotites. 
Overall, most of the SWIR peridotites have been subjected to 
ancient partial melting events that are much older than the spreading 
age of the SWIR, which initiated in the Jurassic (26). Our results 
imply that the oldest model ages identified in the SWIR are domi- 
nantly Archean west of the Andrew Bain transform fault and late 
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Proterozoic in age to the east. Moreover, the Knorr peridotites 
show variable textures and model characteristics (see Materials and 
Methods), and therefore, they cannot have come from a single out- 
crop but must represent sampling consistent with the 3.6-km length 
of the dredge track (fig. S9). Therefore, our results substantiate the 
occurrence of an Archean cratonic mantle domain with a scale larger 
than 3 km at the western segment of SWIR. 

The refractory compositions and old ages of the HFK peridotites 
suggest that they represent the same low-density lithologies that are 
found in Archean cratonic keels. The possibility that they might 
represent pieces of depleted Archean oceanic lithosphere, although 
it cannot be ruled out, is hard to be tested, as there are no preserved 
analogs that are actually known to be derived from Archean oceanic 
mantle. By contrast, samples of Archean subcontinental lithosphere 
are preserved as mantle xenoliths in kimberlites. As already pointed 
out further above, the chemical and isotopic characteristics of the HFK 
peridotites are notably similar to those of the cratonic lithosphere. 


We tested this statistically by carrying out Student’s t test comparing 
the HFK samples with a population of 183 South African cratonic 
(SAC) peridotites and found that the two populations are indistin- 
guishable in both Os isotopes and model ages (see Materials and 
Methods). For this reason, we assume that the most plausible origin 
of these peridotites is that they are derived from the Archean SCLM, 
and this will be the working hypothesis for the remainder of 
this paper. 

This still leaves the questions of how and when these rocks were 
originally dislodged from the SCLM and entered the asthenosphere. 
The dominant mechanisms at play remain enigmatic. Loss of keels in 
the North China and Wyoming cratons has been attributed to subduc- 
tion of the Pacific plate in the Mesozoic (27, 28), in stark contrast 
to the tectonic setting around southern Africa that experienced no 
subduction in its recent geologic past. An unresolved problem 
encountered with the different tectonic scenarios invoked for the 
removal of Archean lithosphere lies in the fact that much of this 
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Fig. 4. Numerical modeling of proposed pathway for the recycling of Archean cratonic mantle. The dynamics of a plume entraining craton material into the asthenosphere 
that is then transported, due to relative buoyancy, all the way to the mid-ocean ridge. Two snapshots at 29.5 Ma (A) and 76 Ma (B) after simulation initiates. (C) Preserved 
area of cratonic mantle lithosphere during plume impingement. Model parameters used are indicated in fig. S10 in addition to exploring additional parameter space. See 


Fig. 5 for an artistic rendition of the mechanism, including subsequent subduction. 
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Fig. 5. Proposed pathway for the recycling of Archean cratonic mantle. (1) The cratonic mantle is disrupted by rifting and/or a mantle plume and is fragmented into 
the asthenosphere, after which it is transported to a mid-ocean ridge. (2) The recycled cratonic mantle is incorporated into the oceanic lithospheric mantle at a mid-ocean 
ridge during ultraslow seafloor spreading where it can be exposed and sampled. (3) Ultimately, subduction returns the fragments of transported Archean cratonic mantle 
back to the convective mantle from which it was originally extracted. The artwork was designed by Shen-Wen Wang@Institute of Planets Inc. 


Archean lithospheric material consists of ultrarefractory, low-density 
peridotites, which are buoyant relative to the asthenospheric mantle 
(11, 29). Consequently, delamination mechanisms invoking density- 
driven lithospheric loss are problematic. It has been proposed that 
cratonic keels can increase density by metasomatic addition of mafic 
components that would generate sufficient garnet to increase the 
bulk density of the lithosphere (30, 31). Such a mechanism can be 
ruled out in the present case, however, because the Archean HFK 
samples have retained their original ultrarefractory (commonly 
with <1% Al,O3; table 1) and thus low-density properties, despite 
traces of metasomatism from the presence of rare Cpx (fig. $4). This 
indicates that the cratonic mantle can be detached even when it is 
buoyant, thus favoring shear flow over densification as a mecha- 
nism for SCLM loss. We suggest that the most likely agent for such 
shear flow is a mantle plume. Similar uncertainties exist regarding 
the timing of lithospheric destruction. Some models of net crustal 
growth contend that the mass of the late Archean continental crust 
reached at least 70% of the present crustal mass (32-34), which, if 
true, would imply that large-scale destruction of Archean crust and 
its associated SCLM must have taken place, given its small fraction 
of preserved surface ages today. Therefore, one might expect to find 
pieces of Archean SCLM in the modern mantle sampled at many 
ocean ridges that contain outcrops of mantle peridotites. However, 
the previous record of oceanic peridotites contains no such exam- 
ples of recycled Archean SCLM, perhaps indicating that mantle 
mixing and homogenization processes may have destroyed most or 
all of this evidence. 

Because neither the timing nor the mechanisms of removal of 
the peridotites from the subcontinental lithosphere can be inde- 
pendently constrained, we propose a model for this process that is 
consistent with the regional geology and seismic tomography. We 
suggest that the HFK samples were dislodged from the base of the 
nearest craton, the Kaapvaal craton ~1500 km away, in the relatively 
recent geologic past. There are at least two ways for the Kaapvaal 
Archean cratonic mantle to have entered the asthenosphere. It 
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could have been disrupted along this margin during the Jurassic 
breakup of Gondwana and then transported in the convecting 
asthenosphere (10). Alternatively, or in combination with rifting, 
sustained plume activity beneath the western Gondwana cratons 
since the Late Cretaceous could have weakened the SCLM through 
warming and/or metasomatism, resulting in its erosion and thinning 
(30). Plume-triggered erosion is also supported by studies of mantle 
xenoliths entrained in the older (>110 Ma) and younger (<100 Ma) 
kimberlites in the southwest Kaapvaal craton, which revealed a 
thinning of ~30 km in thickness of the cratonic lithosphere over the 
short time span separating these two eruptive events (35). Trans- 
porting the eroded Archean lithospheric mantle from the Kaapvaal 
craton to the SWIR also fits well with the current seismically inferred 
mantle convective flow field (Fig. 1), which appears to be dominated 
by a major upwelling originating in the lower mantle beneath 
southern Africa and radiating away from Africa along the base of the 
asthenosphere under the Afar plume (36). Seismic tomography of 
the mantle beneath Africa (37) shows one branch of the sub-African 
“plume tree” extending toward the Bouvet triple junction, and it 
might be this branch, or its precursor, the Karoo plume, which 
eroded the base of the Kaapvaal craton and carried the lithospheric 
fragments toward the spreading ridge. 

We have numerically modeled SCLM delamination by simple 
shear flow resulting from the impingement of a mantle plume on 
the base of the lithosphere and recycling of the dislodged SCLM to 
a spreading ridge (see Fig. 4 and Materials and Methods). The 
dislodgement process is similar to a previous conceptual model for 
lithospheric detachment (38), although our model does not require 
partial melting to occur at the base of the lithosphere. The complete 
recycling process is capable of shearing off sheets of lithospheric 
material, entraining it, introducing it into the asthenosphere, and 
ultimately returning parts of it to a spreading ridge ~1500 km away 
as in our map area (Fig. 1A). The numerical modeling results 
(Fig. 4) thus confirm the feasibility of the recycling mechanism that 
we propose to explain our data (Fig. 5). Our finding that “pristine” 
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ancient lithospheric material appears to exist in at least one portion 
of the global ridge system should reopen the search for other such 
occurrences and encourage a reexamination of the dynamics of 
such a recycling process. The apparent scarcity of the survival of 
such material in the present-day oceanic mantle might be caused by 
efficient convective stirring (4), which may gradually eradicate the 
memory of such ancient processes. Alternatively, it might simply be 
the result of undersampling of oceanic peridotites, which have been 
examined for their possible “Archean memory” by relatively few 
studies so far (3-5, 39). 

If such cratonic recycling can occur today, it has likely occurred 
throughout Earth history. Considerations from a mantle perspective 
of net continental growth (i.e., accounting for the effects of conti- 
nental recycling and reworking) suggest that most or all of the 
present-day continental mass is likely to have been reached by the 
early Archean (32, 34). However, little of this ancient crust and its 
lithospheric keel is preserved today in Earth’s surface ages. We 
propose that much of the cratonic mantle and the Archean crust it 
protects has likely been subducted into the convecting mantle after 
being recycled to mid-ocean ridges. Our modeling indicates that 
this proposed mechanism is highly efficient, removing ~20% of 
the SCLM adjacent to the mantle plume in ~100 million years (Ma) 
(Fig. 4C), where intermittently integrating such a rate of removal 
during repeated episodes of rifting since the early Archean would 
have cumulatively removed a substantial volume of cratonic mantle 
to account for the observed present-day deficiency. Cratonic litho- 
sphere delamination and recycling caused by plume-induced mantle 
shear flow may have been a common mechanism for craton de- 
struction in Archean time, especially if mantle plumes were the 
primary agent for generating widespread, hot Archean komatiites 
(40). The history of crustal growth and destruction has been widely 
studied by geochemical box modeling (34) and by examining the 
surviving parts of the ancient crust (41). We suggest that the key for 
understanding the destruction part of this story is most likely not 
preserved in the surviving crust but must be sought in the mantle. 


MATERIALS AND METHODS 

Materials: SWIR abyssal peridotites 

The SWIR extends 7700 km between the Bouvet and Rodriguez 
triple junctions (Fig. 1A) and is one of the slowest spreading mid- 
ocean ridge with a near-constant full spreading rate of ~14 mm year’. 
Much of the SWIR has a very oblique (~60°) trend from the plate 
spreading direction, reflecting the dominance of numerous large 
transform offsets (42). The SWIR may be divided into a number of 
subsections based on changes in the obliquity of the ridge axis and 
the variation of regional axial depths. Between the Shaka FZ at 9°E 
and the Du Toit FZ at 25°E, with ~1000 km in length, it consists of 
two sharply contrasting supersegments (42), i.e., the orthogonal 
supersegment (16°E to 25°E) and the oblique segment (9°E to 16°E). 
The ridge in the orthogonal segment is oriented approximately 
perpendicularly to the spreading direction and is linked by a 150-km 
slightly oblique (~10°) section. The oblique supersegment begins at 
the 16°E discontinuity with a 500-m drop in rift valley depth. The 
ridge in this supersegment is highly oblique (50° to 60°) and has a 
mean depth of ~4000 m. Three amagmatic accretionary segments 
exist within the oblique supersegment (i.e., 9°30’E to 11°E, 11°35’E 
to 14°15’E, and 14°54’E to 15°45’E), which are connected by the 
Narrowgate segment and the Joseph Mayes seamount. The three 
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amagmatic segments have an average axial depth of 4151 + 425 m 
and variable rift valley widths of ~30 to 50 km (43), from which 
abundant peridotites, scattered basalts, minimal diabases, and 
gabbros have been recovered. The Narrowgate segment (11°45’E to 
14°54’E) is a volcanically robust, orthogonal spreading segment 
with a ~2700-m deep and ~20-km-wide rift valley. The Joseph Mayes 
seamount between 11°E and 11°35’E is spreading approximately 
orthogonally and is perhaps the largest non-hot spot volcano on 
the global ridge system (43). 

At 25°E, the SWIR is abruptly offset about 1800 km to the northeast 
by the four FZs Du Toit, Andrew Bain, Marion, and Prince Edward. 
The three subsections between these four FZs display almost con- 
stant obliquities of 25° and have an average axial depth of ~3200 m 
(44). The Andrew Bain FZ has been suggested as the transition 
between Indian-like and Atlantic-like mantle sources (45). From the 
Andrew Bain FZ to the Gallieni FZ, the SWIR crosses the flank of 
the Southern Ocean geoid high, following a subdued positive arch 
over the Marion swell, starting at a depth of ~4800 m at 27°E before 
rising abruptly to a depth of 858 m at 36°E and then down to 5150 m 
at 63°E. The SWIR east of the Indomed FZ was formed by propaga- 
tion of the Rodriguez triple junction (RTJ) over the past 64 Ma. It 
has regularly spaced nontransform discontinuities, short oblique 
amagmatic segments, and the Atlantis I], Novara, and Melville 
transforms. From 61°E to the RT), like in the 9°E to 25°E section, 
the SWIR is devoid of long-lived transform and nontransform 
discontinuities. 

Our samples come from different regions along the SWIR. One 
suite, dredged by RV Knorr Cruise 162-9 (Dredge 37), is from the 
western end of the ~190-km-long highly oblique amagmatic segment 
(11°38.7’E to 14°17.2’E; Fig. 1B), ~25 km from the large alkalic 
Joseph Mayes seamount (43). The obliquity of the amagmatic 
segment reduces the local mantle upwelling rate, resulting in con- 
ductive cooling capping the melting column at great depth limiting 
melting, likely accounting for the enriched volcanism to the west 
(42). Dredge 162-9 was made from 3907 to 3433 m from near the 
base of the northern wall of the 9°E to 16°E rift valley (fig. S9) along 
a 3.6-km-long track over rough terrain. The dredge recovered a 
diverse assemblage mainly consisting of partially serpentinized and 
weathered peridotites (249 kg), with minor 0.4-kg gabbro (0.4 kg), 
diabase (4 kg), basalt (1 kg), and volcaniclastics (1.5 kg). The 
peridotites show variations in petrographic textures (fig. $10) and 
degrees of alteration (fig. S11). Another suite was dredged ~1700 km 
to the east across the Andrew Bain FZ by RV Melville Protea 
Expedition Leg 5 from the Marion and Prince Edward FZs, respec- 
tively (Fig. 1C), Le., PROT 5-19 (46°34’S, 33°49’E) and PROT 5-24 
(45°33’S, 35°10’E). 


Major and trace elements 

All measurements reported here were conducted at the Institute of 
Geology and Geophysics, Chinese Academy of Sciences (IGGCAS). 
Whole-rock major elements were analyzed using x-ray fluorescence. 
Mineral major elements were measured by a Cameca SX Five 
electron microprobe, using an acceleration potential of 15 kV, a 
beam current of 30 nA, and a spot size of 5 um. Trace elements of 
both Cpx and Opx were determined by laser ablation system, which 
consists of an Element XR HR-ICP-MS instrument (Thermo Fisher 
Scientific, USA) and a 193-nm ArF excimer Analyte G2 laser system 
(Teledyne CETAC Technologies, USA). The isotopes were measured 
using a peak-hopping mode with a beam size of 60 1m and a repetition 
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rate of 5 Hz. *’Si was used as an internal standard. NIST SRM 612 
reference glass was used for external calibration. Both BCR-2 and 
BIR-1 glass were used for quality control monitoring. 


HSEs and Re-Os isotopes 

Whole-rock Re-Os isotopes were analyzed at the IGGCAS using the 
isotope dilution method (46). About 2 g of powders, together 
with a '*’Re-'°Os spike and a mixed '?'Ir-’Ru-!*Pt-'Pd spike, 
was mixed in Carius tubes, which were digested with reverse aqua 
regia (i.e., 3 ml of 12 N HCl and 6 ml of 16 N HNOs) at 240°C for 
~72 hours in an oven. Osmium was extracted from the aqua regia 
solution by solvent extraction into CCl, and further purified by 
microdistillation (47). Other HSEs were first separated from the 
solution into subgroups (Re-Ru, Ir-Pt, and Pd) using a 2-ml anion 
exchange resin (AG-1 x 8, 100 to 200 mesh). Subsequently, the 
Re-Ru was further purified by a 0.25-ml anion exchange resin; the 
Pd and Ir-Pt were further purified by an Eichrom LN spec resin to 
completely remove Zr and Hf. 

Osmium concentrations and isotopes were measured by negative 
thermal ionization mass spectrometry on a GV IsoProbe-T instru- 
ment in static mode using Faraday cups. A Ba(OH), solution was 
used as an ion emitter. The measured Os isotopes were corrected for 
mass fractionation using the '”’Os/'**Os ratio of 3.0827. The Nier oxygen 
isotope composition ('70/'SO = 0.0003708 and '80/'°O = 0.002045) 
has been used for oxide correction. The in-run precisions for Os 
isotopic measurements were better than 0.2% (20) for all the 
samples. The Johnson-Matthey standard of UMD was used as an 
external standard, yielding an '*’Os/'**Os ratio of 0.11378 + 2 (20; 
n = 5). Concentrations of other HSE were measured on a Thermo 
Fisher Scientific Neptune MC-ICPMS with an electron multiplier 
in peak-jumping mode or using Faraday cups in static mode, de- 
pending on the measured signal intensity. In-run precisions for 

*Re/*/Re, 'r/! Ir, PRu/™'Ru, 4Pt/!°°Pt, and '°°Pd/'°°Pd were 
0.1 to 0.3% (26). The Re, Ir, Ru, Pt, and Pd standards were used 
to correct mass fractionations. The total procedural blanks were 
3 + 1 pg for Os, 5 + 2 pg for Re, 1 + 0.5 pg for Ir, 10 + 5 pg for Ru, 
20 + 10 pg for Pt, and 10 + 5 pg for Pd (n = 10). An '8’Os/'*"Os ratio 
of ~0.15 was obtained for the blank. The standard UB-N was used 
to monitor the accuracy of the analytical procedure; results are 
given in table S4. 


Age estimate of mantle peridotites using Re-Os isotopes 

Multiple methods have been proposed to use the Re-Os isotope 
system to determine the ages of peridotites (22, 48). As for most 
chronometers, the isochron is expected to provide a robust age 
estimate. However, meaningful Re-Os isochrons have rarely been 
obtained for mantle peridotites, probably because of the disturbance 
of this system by late-stage processes after melting in mantle 
peridotites (e.g., melt metasomatism). Lacking reliable isochrons, 
two different kinds of Os model ages have been widely applied for 
mantle peridoites (49). One is the Os model age, Ta, which is 
calculated by using the measured '*’Re/'**Os ratio of the sample to 
extrapolate the measured '*’Os/'*8Os ratio back in time until it 
intersects the evolution curve of a model mantle reservoir, e.g., the 
chondritic upper mantle (49) or the PUM (20). However, the 
reliability of the Tma age is often compromised by Re mobility in 
mantle peridotites (the same problem encountered in determining 
Re-Os isochrons). For this reason, the most widely used Os model 
age is the rhenium depletion model age, Trp, which is calculated by 
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directly comparing the measured '*’Os/1*Os ratio of the sample to 
the evolution curve of a model mantle reservoir (49). As the Trp age 
ignores the in-growth of '*’Os due to Re in mantle residues after 
melting, it provides a minimum estimate of the true melt-depletion 
age. Therefore, the Trp age is most applicable for strongly refractory 
mantle peridotites in which Re is almost completely removed 
during high degrees of partial melting, with little Re remaining in 
the mantle residues. It should be noted that Trp ages calculated 
relative to the evolution curve of the PUM model mantle [current 
1879)5/188O5 = 0.1296; (20)], as shown in table $4, are ~200 Ma 
older than those calculated relative to “chondritic” mantle [current 
1875/1885 = 0.127; (49)]; however, Archean ages are found for 
some Knorr samples regardless of which evolution curve is used. 


Thermomechanical numerical modeling 

The two-dimensional thermomechanical numerical code I2VIS 
based on the finite difference method with marker-in-cell technique 
(50) was used to solve the mass, momentum, and energy conserva- 
tion (Eqs. 1 to 3) using the Boussinesq approximation 


V-V=0 (1) 
-VP+V[n(VV + VV')]=pg (2) 
pCy or = V-(kVT)+H, +H, + Hy (3) 


where V is the velocity, P is the dynamic pressure, p is the density, 
g is the gravity acceleration, C, is the specific heat, T is the tem- 
perature, k is the thermal conductivity, H, is radioactive heating, 
H, = ty x €y is shear heating, and H, = Tapgvy is adiabatic heating. 
The effective viscosity is given by y = (1/Naisa + 1/naite) 2, where nist 
and ngif¢ are dislocation and diffusion creep viscosities, respectively. 

The effective density of rocks depending on the temperature and 
pressure is given by 

P = po(1 — a(T — To)) (1 + BP — Po)) (4) 
where Po is the density of rock at the reference temperature 
(To=273 K) and pressure (Po=0 Pa), o is thermal expansivity, and 
B is compressibility. 

The model domain is 2700 km by 675 km in x-y directions, 
which is resolved by 901 x 226 nodes and more than 7 million 
markers to trace physical properties (a relatively higher resolution 
with 3 km by 1.5 km is also tested, and invisible difference is ob- 
served; fig. $12). All boundaries are using free slip condition, and to 
mimic a free surface in the finite difference, we impose a 15-km-thick 
sticky air (10!8 Pa-s) (51, 52). The lithosphere has a linear tempera- 
ture profile with 0°C at the surface and 1350°C at the bottom, while 
the asthenospheric mantle has an adiabatic thermal gradient of 0.5 K 
km™'. The top and bottom boundaries are isotherm to have fixed 
temperatures of 0°C at the surface and 1680°C at the bottom, re- 
spectively. The side boundaries have zero heat flux. 

The laboratory-derived parameters of the flow law for olivine 
(53) are listed in table S5. The initial model setup is composed of a 
200-km-thick cratonic lithosphere and a 100-km-thick orogenic 
lithosphere, and the mid-ocean ridge is located 1000 to 2500 km 
away from the craton. The thermal age of the oceanic lithosphere 
linearly increases from 0 Ma at the ridge to 30 Ma at 200 km away, 
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while a half-space cooling model is used for its thermal structure. 
We assume that a hot mantle plume consistently sits below the craton. 
Its radius and excess temperature are varied in our model (fig. $13) 
to investigate their influences on the ability of entrained cratonic 
material to be transported to the far-field ridge. Because of its low 
thermal structure, the cratonic lithosphere has a large effective viscosity 
(>107° Pa:s) to resist persistent mantle convection. However, mantle 
upwelling in a strong and large plume can erode the lithosphere at 
the craton’s edge where thermal structures are remarkably different 
between two terranes and edge-driven convection is pervasive (54). 
The arrival of the mantle plume to the base of the lithosphere tends 
to be tilted to the orogenic lithosphere side (55), and in conse- 
quence, the continuously fed mantle upwelling entrained the eroded 
or delaminated cratonic mantle to a distance of more than 1500 km 
to reach the ridge in only a few tens of million years to 100 Ma. 


Student's t test 

Student’s t test was conducted to determine whether the HFK 
peridotites and the SAC peridotites represent similar or distinct age 
populations. Because there are more SAC (n = 183) than HFK 
(n = 20) peridotites and Student’s t test must be conducted on two 
populations of equal sample size, we resampled 20 ages from the 
SAC peridotites 100 times. Given the skewed age distributions of 
both datasets (Fig. 3), the t test was conducted with the one-tailed 
distribution. Because the two sample sets are independent of each 
other, the two-sample f test (assuming equal variance) was used. 
The average P value of the 100 t tests from the 100 resamplings was 
calculated. Far exceeding the typical threshold of 0.05, the value of 
P = 0.26 indicates that the null hypothesis that the two age popula- 
tions are the same cannot be rejected. This result is consistent with 
the HFK peridotites being sampled from the same age population 
of the SAC peridotites. If the HFK peridotites were not sampled from 
the SAC peridotites, their statistically identical age populations would 
have to be purely coincidental. 


SUPPLEMENTARY MATERIALS 


Supplementary material for this article is available at https://science.org/doi/10.1126/ 
sciadv.abn6749 
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